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EQUATIONS  OF  FLOW  IN  A  RAREFIED  ATMOSPHERE* 


Harold  Grad 
ABSTRACT 

The  problem  of  the  aerodynamic  or  flow  properties  of  the 
upper  atmosphere  is  characterized  by  the  presence  of  a  large 
number  of  simultaneous  and  interacting  phenomena  and  a  vast  range 
of  values  taken  by  the  relevant  physical  parameters  with  an  en- 
suing bewildering  array  of  different  flow  regimes.   As  a  working 
rule,  one  can  hope  to  develop  a  non-trivial  mathematical  theory 
if  at  most  one  physical  parameter  (e.g.  the  Mach  number)  is 
allowed  to  range  over  all  its  values,  and  all  other  parameters 
(e.g.  the  mean-free-path)  are  taken  to  be  either  very  small  or 
very  large.   In  the  first  part  of  this  paper,  we  present  a  survey 
of  some  of  the  more  interesting  limiting  cases,  and  we  mention  only 
briefly  the  multitude  of  transition  regimes.   In  the  second  part  we 
consider,  as  an  illustrative  problem,  the  "free  flow"  around  an 
object  much  smaller  than  the  mean-free-path.   Surprisingly, 
whether  or  not  this  problem  should  be  handled  by  the  traditional 
free  molecular  flow  techniques,  or  by  continuum  techniques,  or  as 
a  problem  in  transition  flow  turns  out  to  depend  on  what  question 
has  been  asked.   In  particular,  one  can  find  in  this  single 


* 
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problem  all  the  features  of  the  various  flow  regimes  for  arbitrary 
Knudsen  number  as  well  as  Mach  number.   The  asymptotic  flow  at 
infinity  can  be  described  by  an  appropriate  solution  of  the  macro- 
scopic non-dissipative  Suler  equations  including  a  wake  behind  the 
body  and  a  shock  wave  in  the  supersonic  case.   The  fine  structure 
of  the  wake  is  similar  to  the  classical  Oseen  solution.   Exactly 
the  same  asymptotic  behavior  is  obtained  for  the  flow  around  a 
body  of  arbitrary  size  (not  in  free  flow).   Thus  the  asymptotic 
flow  depends  on  Mach  nurriber  only  and  not  on  Knudsen  number. 

1.   SURVEY  OF  THE  FLOW  REGIMES 

In  the  problem  of  the  flow  of  an  ordinary  (unionized)  gas 
around  an  object  there  are  essentially  two  relevant  parameters, 
the  mean-free-path  (Knudsen  number)  and  the  speed  (Mach  number). 
The  two  extremes  of  small  and  large  m.ean-free-path  have  been  ana- 
lyzed in  some  detail.   In  the  former  case  we  choose  one  of  the 
conventional  formulations  of  continuum  gas  dynamics,  the  sim- 
plest being  the  non-dissipative  Euler  equations.   For  the  latter 
we  solve  the  simple  kinetic  theory  problem  of  non-interacting 
molecules  which  follow  straight  line  paths.   The  only  complica- 
tion is  the  selection  of  an  appropriate  boundary  condition. 
Varying  the  flow  speed  presents  no  difficulty  in  the  free-flow 
extreme.   It  presents  considerable  difficulty  in  the  continuxm 
formulation,  but  there  exists  an  extensive  literature  which 
covers  this  subject  in  detail.   The  transition  regime  with  re- 
gard to  mean-free-path  is  much  less  understood  than  the  transl- 
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tion  with  regard  to  Mach  number  and,  for  the  flow  around  a  finite 

object,  has  hardly  been  treated  except  in  the  low  speed  limit 

(see  Refs.  (l),(2);  also  (5),(^)  and  additional  Refs.  quoted  there). 

The  introduction  of  ionization  opens  a  veritable  Pandora's 
box  of  complex  phenomena.   To  start,  consider  the  mean-free-path 
which  is  no  longer  unique  and  does  not,  by  itself,  give  a  criterion 
for  either  free  flow  or  simple  continuum  flow.   There  is  a  mean- 
free-path  for  neutral-neutral  encounters,  for  encounters  between 
charged  particles,  and  for  encounters  between  a  neutral  and  a 
charged  particle.   For  simplicity  let  us  ignore  inelastic  effects 
such  as  ionization,  recombination,  charge  exchange,  and  related 
surface  effects.   We  speak  of  free  paths  rather  than  collision 
times  since  the  former  are  about  the  same  size  for  electrons  and 
for  ions.   We  also  choose  to  ignore  as  minor  the  difference 
between  neutral -neutral  and  neutral -charged  cross-sections.   This 
leaves  two  distinct  mean-free-paths,  the  neutral  and  the  Coulomb, 
At  the  temperatures  existing  in  the  atmosphere,  the  Coulomb  cross- 
section  is  much  larger  than  the  neutral  cross-section.   This  is 
counterbalanced  by  the  fact  that  the  neutral  density  is  usually 
much  larger  than  the  charged  particle  density.   The  Coulomb  mean- 
free-path  is  always  large  (at  least  300  meters),  but  the  neutral 
mean-free-path  can  range  from  microscopic  values  at  low  altitudes 
to  values  large  compared  to  the  Coulomb  mean-free-path. 

In  contrast  to  the  situation  for  neutral  gas,  we  shall  see 
that  neither  does  a  small  neutral  mean-free-path  (concomitant  with 
a  small  mean-free-path  for  charged-neutral  encounters)  validate  a 
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simple  continuum  theory,  nor  do  large  values  of  both  the  neutral 
and  Coulomb  mean-free-paths  imply  that  free  flow  prevails.   The 
explanation  lies  in  the  long-range  nature  of  Coulomb  forces.   In 
one  aspect,  the  relevant  dimension  is  the  Debye  shielding  length, 
d.   For  our  purposes,  the  significance  of  this  length  is  that 
removal  of  all  the  charged  particles  of  one  sign  within  a  cube 
of  side  d  will  create  an  appreciable  electric  field  (i.e.  a 

potential  comparable  to  kT).   The  charge  removed  is  nd  e  and  the 

3  2 
potential  at  a  distance  d  is  nd  e  /x:  d  which,  set  equal  to  kT, 

yields  the  result 

(1.1)  d^  =  ^  . 

ne 

If  the  object  is  much  smaller  than  the  Debye  length  as  well  as 
the  mean-free-path,  then  the  perturbation  of  the  flow  induced  by 
the  object  includes  a  negligible  electric  field  perturbation,  and 
the  charged  particle  orbits  can  be  computed  as  in  free  flow.   On 
the  other  hand,  if  the  Debye  length  is  small  compared  to  macro- 
scopic dimensions,  then  we  can  conclude  that  the  ion  and  electron 
densities  are  almost  equal  everywhere  in  the  flow,  and  we  can  use 
greatly  simplified  single-fluid  macroscopic  equations.   Of  course, 
this  quasi-neutrality  of  the  medium  does  not  imply  absence  of 
electric  fields;  on  the  contrary,  any  but  enormous  fields  can 
arise  by  only  very  slight  deviations  from  neutrality. 

In  the  transition  regime  (finite  Debye  length)  we  have  a 
so-called  self-consistent  field  problem.   The  particle  orbits  are 
to  be  computed  in  an  assumed  electric  field  which  must  them  be  made 
consistent  with  the  solution  of  Poisson's  equation  in  which  the 
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charge  source  is  determined  by  the  particle  motions   (e.g.  see 
Ref.  (5)).   For  a  typical  density  of  10^  electrons  per  cm  in 
the  ionosphere,  the  Debye  length  is  on  the  order  of  several  milli- 
meters. 

An  analogous  effect  arises  from  the  magnetic  field.   Taking 
the  maximum  current  density  to  be  that  resulting  from  all  particles 


of  one  sign  moving  at  thermal  speed,  J  =  ne  /kT/m  ,  we  compute 
a  mlnimiom  transition  distance,  6,,  over  which  an  appreciable  mag- 
netic field  can  be  produced  by  setting  the  resulting  magnetic 

energy  B  /2ii  equal  to  the  gas  energy  p  where  J -^  B/p.6-,.   We 

2  '' 

obtain  (using  M.?t=  l/c  ) 

(1.2)  6/  =  lB£%'=^d2  . 

ne 

If  6^  is  large  compared  to  the  size  of  the  object,  then  the  mag- 
netic field  perturbation  can  be  ignored  and  we  have  free  flow 
relative  to  the  magnetic  field.   If  5-j^  is  small  compared  to  macro- 
scopic dimensions,  then  we  conclude  that  the  mean  flow  speeds  of 
the  ion  and  electron  gases  are  approximately  equal.   Even  very 
large  currents  are,  in  this  case,  easily  produced  by  small  relative 
velocities.   In  the  transition  regime,  one  must  compute  the  flow 
of  the  two  fluids  in  an  assvuned  magnetic  field  which  is  to  be  con- 
sistent with  the  solution  of  Maxwell's  equations  containing  the 
related  current  source.   From  the  formula  (1.2),  it  is  clear  that 
the  electric  field  criterion  for  free  flow  is  the  more  stringent, 
whereas  the  magnetic  criterion  for  continuum  flow  is  the  more 
stringent.   A  typical  value  of  6,  in  the  ionosphere  is  about  10 
meters. 

-  7  - 


In  the  presence  of  an  external  magnetic  field,  B  ,  one  can 
compute  an  alternative  distance,  6p,  as  the  mlnlmiim  distance  over 
which  the  Induced  magnetic  field  can  be  comparable  to  the  external 
field.  We  obtain 

(1.3)  5^  =  5^/  /T  =  a/3 

where 

is  a  mean  Larmor  radius  and 
(1.5)         p  -  — 2- — 

is  the  ratio  of  gas  to  magnetic  energy.   In  the  presence  of  an 
external  field,  5  can  also  be  compared  to  the  Larmor  radius, 
5,  =  A/v/p~.   In  the  special  case  of  equipartltlon  of  magnetic 
and  kinetic  energy,  p  =  1,  we  have  5  =  6p  «  A. 

In  general,  we  distinguish  6,  ,  which  is  the  minlmijim  distance 
required  for  Induced  currents  to  have  an  effect  on  the  gas,  from 
5p,  which  is  the  distance  over  which  there  can  be  an  effect  on 
the  ambient  magnetic  field.   In  the  ionosphere,  the  energy  density 
of  the  earth's  field  is  about  10  times  the  kinetic  energy  of  the 
ionized  component,  and  it  becomes  comparable  to  the  total  gas 
energy  at  about  135  1<to« 

There  is  another  case  similar  to  the  ones  described  by  the 
lengths  5-,  and  5-.   Suppose  p  is  small  and  6,  is  large  (as  in  the 
ionosphere).   The  induced  currents  are  not  large  enough  to  perturb 
the  ambient  magnetic  field,  nor  can  they  create  a  large  enough 
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magnetic  pressure  to  be  comparable  to  the  gas  pressure,  p. 
However,  the  force  on  the  ionized  component  Is  measured  by  the 
product  of  the  small  Induced  current  and  the  relatively  large 
ambient  field,  B  .   The  criterion  that  this  force  be  appreciable 
when  compared  to  p  yields  a  length  6^  =  ^.   The  mean  electron 
Larmor  radius  can  have  a  value  from  one  to  several  centimeters 
In  the  Ionosphere. 

The  relative  sizes  of  the  electric  and  magnetic  Lorentz 
forces,  qE  and  J  x  B  ,  depend  on  whether  the  electric  field  Is 
governed  by  the  Debye  length  (as  In  the  boundary  layer  surround- 
ing the  body)  or  by  the  size  of  the  object  (in  the  main  flow). 
In  the  boundary  layer,  the  electric  force  dominates  by  the  ratio 
A/d  ,  whereas  in  the  main  flow,  the  magnetic  force  dominates  by 
the  ratio  of  the  cube  of  a  macroscopic  dimension  to  the  product 

The  next  length  we  consider  is  the  Larmor  radius  In  Its  prim- 
itive interpretation  in  terms  of  particle  orbits.   In  the  free 
flow  regime,  this  offers  complications  very  similar  to  those 
existing  in  conventional  free  flow  problems  in  which  there  occur 
multiple  reflections  with  the  boundary.   One  obtains  an  Integral 
equation  for  the  incident  mass  flow  on  the  surface  of  the  object. 

In  the  continuum  regime,  the  comparative  sizes  of  Larmor 
radius  and  mean-free-path  (both  small)  enter  into  the  size  of 
the  Hall  effect  term  in  Ohm's  law.   Ohm's  law  can  be  written  in 
the  form 
(1.6)  E  +  uxB  =  R|^J,,  +i^JxB  +  Rj.Jl 
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where 


(1.7)  J(/  =  -\   (J-B)B  ,   Jo.  =  4r  (B  X  J)  X  B. 

B  B 


J  =  J(i  +  Ji. 

The  resistivity  is  essentially  proportional  to  the  collision 
frequency  or  reciprocal  mean-free-path,  while  the  Hall  coeffi- 
cient, >|_B,  is  proportional  to  the  Larmor  frequency  or  reciprocal 
Larmor  radius.   For  a  fully  ionized  two -component  gas  (one 
species  of  ion)  the  Hall  effect  arises  as  a  consequence  of  choosing 
the  ''total"  fluid  velocity,  u,  as  the  mass  flow  velocity;  the 
effect  would  not  be  present  (at  least  not  in  the  expression  of 
Ohm's  law)  if  it  were  taken  as  the  arithmetic  mean  of  the  ion 
and  electron  gas  speeds. 

The  parameter  R,,  is  independent  of  B  while  Rj.  and  t^  depend 
only  slightly  on  B;  (they  may  vary  on  the  order  of  a  factor  two 
when  B  varies  from  zero  to  infinity).  For  a  fully  ionized  two- 
component  gas  Rj,  is  found  to  be  equal  to  R,,  when  thermal  diffusion 
(i.e.  thermoelectricity)  is  ignored.  Inclusion  of  thermal  dif- 
fusion or  a  mixture  of  ionic  species  yields  unequal  values  for 
R,i  and  Rj_. 

There  has  been  a  certain  oversimplification  in  the  statement 
that  continuum  flow  (by  which  we  mean  a  conventional  one-fluid  set 
of  equations  combined  with  Maxwell's  equations  and  Ohm's  law) 
results  when  the  mean-free-paths  together  with  d  and  6  can  be 
considered  to  be  small.   One  usually  considers  the  conditions 
that  mean-free-path,  Larmor  radius,  and  Debye  length  be  macro- 


scopically  small  to  be  roughly  equivalent  to  the  requirements 
that  collision  frequency,  Larmor  frequency,  and  plasma  frequency 
respectively  be  large  compared  to  relevant  macroscopic  frequencies, 
Implied  in  this  equivalence  is  the  existence  of  a  macroscopic 
speed  which  relates  a  length  to  the  appropriate  frequency.   This 
is  not  so  for  the  case  of  plasma  oscillations  which  have  the 
peculiar  property  that  the  frequency  rather  than  the  phase  speed 
is  approximately  independent  of  wavelength.   In  particular,  it  is 
possible  to  consider  plasma  oscillations  as  being  in  a  sense 
macroscopic,  i.e.  with  a  macroscopic  wave  length  much  larger  than 
the  Debye  length,  even  though  the  frequency  is  usually  too  high 
to  be  considered  within  the  province  of  fluid  dynamics.   Plasma 
oscillations  are  incorporated  into  a  macroscopic  theory  in  which 
the  Debye  length  has  been  removed  by  taking  an  inertial  form  of 
Ohm's  law  in  which  the  time  derivative  bJ/dt  is  kept. 

Summarizing,  we  would  recommend  that  a  macroscopic  single- 
fluid  theory  be  used  to  describe  the  flow  near  an  object  when 
the  mean-free-paths  and  lengths  d  and  5  are  small  compared  to 
the  size  of  the  object.   There  is  a  range  (up  to  100  km)  where 
the  neutral  mean-free-path  is  small  compared  to  a  body  size  of, 
say,  one  meter.   The  charged  particle  flow  would  be  dominated  by 
the  charged -neutral  mean-free-path  which  is  about  the  same  size 
as  the  neutral  mean-free-path.   The  Debye  length  and  the  electron 
Larmor  radius  (6,)  are  small,  but  6,  is   fairly  large  and  6„  is 
very  large.   One  can  therefore  use  a  simplified  version  of  the 


continuum  equations  In  which  the  earth's  magnetic  field  Is  con- 
sidered to  be  given  and  Is  unaffected  by  the  object,  but  the 
magnetic  force  on  the  charged  component  cannot  be  Ignored.   However, 
since  the  neutral  gas  pressure  Is  much  larger  than  the  Ionized 
gas  pressure,  this  magnetic  force  can  be  neglected  In  the  total 
flow  equations  and  enters  only  as  an  additional  Hall  term  in 
Ohm's  law. 

Above  some  I50  km,  all  mean-free-paths  are  large.   Since  the 
Debye  length  is  still  small,  one  must  use  a  self-consistent  elec- 
tric field  analysis  for  the  flow  of  the  charged  components.   The 
magnetic  forces  can  be  important  in  perturbing  the  flow,  but  they 

arise  in  the  form.  J  x  B  where  B  is  the  unoerturbed  magnetic 

o         o  ' 

field.   From  our  assumptions  we  would  conclude  that,  in  this 
region,  the  neutral  particle  flow  and  the  charged  particle  flow 
could  be  solved  Independently.   However,  this  is  true  only  Insofar 
as  they  are  not  related  by  surface  effects  and  by  inelastic 
collisions. 

2.   FREE  FLOW  REEXAMINED 

Consider  the  free  flow  of  a  neutral  gas  about  a  sphere.   The 
conventional  solution  can  be  described  as  follows.   Given  a  point 
In  space,  x,  and  a  molecular  velocity,  ^.   The  number  density  of 
particles  (i.e.  the  distribution  function)  is  the  same  as  at 
infinity  if  a  line  traced  backwards  with  this  velocity  does  not 
Intersect  the  sphere;  the  n\imber  density  is  equal  to  the  number 
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of  particles  leaving  the  sphere  with  this  velocity  if  the  path 
traced  backwards  does  intersect  the  sphere.   Thus  the  distribution 
function  can  be  considered  to  be  the  sum  of  two  distribution  func- 
tions.  First  there  is  the  unperturbed  distribution  function  at 
infinity.   Then  there  is  a  function  which  is  zero  except  on  the 
cone  subtended  by  the  sphere.   This  is  the  difference  between  the 
particles  emitted  by  the  sphere  and  those  shielded  by  the  sphere 
from  the  initial  flow. 

This  description  of  the  flow  implies  that  the  mean-free-path 
has  been  set  equal  to  infinity.   To  be  more  precise,  let  us  consider 
a  case  in  which  the  ratio  a/L  of  the  radius  of  the  sphere  to  the 
mean-free-path  is  very  small  but  not  zero.   The  elementary  free 
flow  solution  described  above  is  approximately  valid  for  a  distance 
r  from  the  sphere  which  is  small  compared  to  L;  this  may  be  many 
sphere  radii.   At  a  distance  comparable  to  L,  the  emerging  stream 
will  begin  to  encounter  the  incident  flow  and  both  will  be  altered. 
The  emerging  stream  will  be  deflected  back  and  eventually  Join 
the  main  stream  while  shielding  the  body  from  the  incident  stream 
which  is  somewhat  deflected  around  the  body.   Clearly,  the  elemen- 
tary free  flow  solution  is  not  uniformly  valid  for  large  distances 

from  the  sphere.   However,  since  the  free  flow  perturbation  decays 

2  2 

as  (a/r)   and  is  small  on  the  order  of  (a/L)   by  the  time  collisions 

become  effective,  one  can  say  that  the  free  flow  solution  is  accu- 
rate in  the  entire  region  where  there  is  a  significant  perturbation 
of  the  original  flow.   Nevertheless,  there  is  some  interest  in 
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determining  the  ultimate  behavior  of  the  flow,  far  from  the 
object.   This  requires  a  full  treatment  of  the  collision  problem 
and  the  Boltzmann  equation.   One  can  expect  free  flow  to  prevail 
near  the  sphere,  continuum  flow  at  a  large  distance,  and  a  complex 
transition  in  between. 

This  question  is  related  to  the  more  complex  problem  of  the 
transition  of  the  flow  near  the  body  from  free  flow  to  continuum 
flow  as  (a/L)  increases.   In  particular,  for  this  problem  one 
would  like  to  see  how  a  detached  shock  (for  supersonic  flow),  and 
a  viscous  boundary  layer,  and  possibly  a  wake  behind  the  body 
develop.   A  striking  result  that  we  find  is  that  the  "detached" 
shock  and  wake  exist  even  in  the  free  flow  limit  of  very  small 
(a/L)  provided  that  we  look  for  them  far  from  the  body.   This  is 
in  contrast  to  the  near  flow  which  does  not  distinguish  at  all 
between  subsonic  and  supersonic  speeds.   Succinctly,  the  near 
flow  sees  no  important  qualitative  dependence  on  Mach  number,  while 
the  far  flow  (we  shall  verify)  sees  no  dependence  at  all  on  the 
Knudsen  n\;miber.   The  simplicity  of  the  asymptotic  behavior  at 
large  distances  from  a   sphere  in  free  flow  is  a  consequence  of 
the  fact  that  the  flow  is  to  a  very  good  approximation  that  de- 
scribed by  the  inviscld  continuum  equations  with  a  point  singularity 
replacing  the  sphere. 

First  let  us  examine  the  asymptotic  behavior  macroscopically. 
To  be  precise,  we  consider  a  scale  of  length  on  which  not  only  a, 
but  even  L  is  small.   There  is  a  general  result  for  the  Boltzmann 
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equation  that  If  we  allow  the  mean-free-path  to  approach  zero  in 
a  given  problem,  the  solution  to  the  Boltzmann  equation  becomes 
locally  Maxwellian  and  the  flow  parameters  satisfy  the  conventional 
non-dissipatlve  equations  of  gas  dynamics  locally,  i.e.  except 
in  certain  narrow  regions  such  as  shock  waves  and  boundary  layers 
(Ref.  (6)).   Of  course,  this  could  be  trivially  true  in  our  problem, 
with  the  appropriate  solution  of  the  fluid  equations  as  the  unper- 
turbed constant  flow.   This  would  be  the  case,  for  example,  if  the 
only  effect  of  collisions  on  the  emerging  stream  were  to  retard 
it  and  produce  a  flow  perturbation  which  decayed,  let  us  say, 
exponentially  rather  than  as  l/r  .   This,  however,  cannot  be. 
However  else  the  sphere  perturbs  the  flow,  it  must  behave  like  a 
momentum  sink;  the  drag  is  a  given  constant  Independent  of  the 
distance  at  which  one  observes.   But  this  singularity  is  macro- 
scopic and  must  make  an  appearance  in  the  momentum  equation.   In 
particular,  we  discover  that  the  macroscopic  solution  decays  at  a 
rate  slower  than  l/r  in  certain  regions. 

On  the  microscopic  level  it  is  possible  to  improve  the  free 
flow  solution  by  taking  into  account  successive  collisions  between 
the  incident  and  emergent  streams,  (see  Ref.  (2),  also  {3),    (^) 
and  additional  Refs.  quoted  there).   All  such  procedures,  whether 
expressed  in  terms  of  iterations  or  as  power  series  in  (a/L), 
seem  to  give  suitable  approximations  only  for  distances  not  too 
far  from  the  body.   This  is  consistent  with  the  actual  picture 
(which  we  shall  verify)  that  the  effect  of  the  emergent  stream 
is  felt  at  a  considerable  distance  away  from  the  body.   Qualita- 
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tively,  one  can  visualize  a  "pile-up"  on  the  order  of  one  mean- 
free-path  In  front  of  the  body  where  the  emergent  stream 
begins  to  deflect  the  Incident  stream  and  Is  Itself  stopped 
and  eventually  turned  back  to  flow  around  the  body.   A  number  of 
mean-free-paths  beyond  the  body,  this  will  have  developed  into  a 
weak  shock  (strength  on  the  order  of  a  /L^ )    if  the  flow  is  super- 
sonic.  From  the  microscopic  picture  it  is  not  at  all  clear  why 
there  should  be  a  radical  change  in  the  behavior  of  the  flow  far 
from  the  object  when  the  speed  changes  from  subsonic  to  supersonic. 
This  distinction  can  only  be  made  some  distance  from  the  body 
after  the  continuum  equations  take  over.   Simiilarly,  directly 
behind  the  body  there  is  a  wake  which  is  only  partially  filled 
in  (residue  on  the  order  of  a^/L")  several  mean-free-paths  down- 
stream, and  its  subsequent  behavior  is  macroscopic. 

In  the  remainder  of  the  paper  we  consider  first  the  elemen- 
tary problem  of  the  linearized  perturbation  induced  by  a  point 
source  using  the  compressible  inviscid   Euler  equations  and  then 
show  that  this  solution  is  the  correct  limiting  flow  far  from 
a  small  object  in  free  flow  obtained  from  the  Boltzmann  equation. 
Of  course,  the  Boltzmann  solution  yields  additional  information 
about  the  flow. 

The  problem  of  the  flow  around  an  object  has  been  studied  in 
great  detail  using  the  Navier-Stokes  equations  (e.g.  Refs.  (?),  (8), 
(9),  (10)).   The  first  solution  we  obtain  (which  is  essentially  the 
fundamental  solution  of  the  linearized  Euler  equations)  is  more 
primitive  and  can  be  obtained  from  the  viscous  solutions  in 
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the  limit  as  the  viscosity  approaches  zero.   The  Boltzmann 
solution  is  more  general  and  can  be  used  to  justify  the  asymptotic 
validity  of  the  non-dissipative  solution  and  also  to  spell  out  the 
exact  regions  of  validity  of  the  Navier-Stokes  solution. 

The  solution  which  we  obtain  for  the  flow  far  from  a  small 
body  in  free  flow  can  be  applied  to  the  asymptotic  behavior  far 
from  a  body  of  arbitrary  size.   This  solution  is,  to  a  certain 
order.  Independent  of  the  Knudsen  number, and  it  is  uniformly 
valid  for  all  Mach  numbers.   As  a  matter  of  fact,  the  asymptotic 
flow  is  completely  determined  by  a  single  parameter,  the  drag  on 
the  body. 

If  ionization  and  electromagnetic  phenomena  are  taken  into 
account,  the  picture  is  much  more  complicated  in  detail,  but  the 
overall  behavior  should  be  quite  similar.   Suppose  we  have  a  body 
which  is  small  compared  to  all  the  relevant  lengths,  viz.  mean- 
free-paths,  Debye  length,  Larmor  radius,  and  magnetic  shielding 
lengths.   The  solution  will  be  an  elementary  free  flow  near  the 
body.   One  passes  through  all  sorts  of  complex  transitional  zones 
corresponding  to  the  reference  lengths  itemized  above.   Far  away 
one  recovers  the  continuum  equations;  moreover,  at  sufficiently 
large  distances  these  equations  will  be  non-dissipative  (including 
infinite  electrical  conductivity).   This  asymptotic  solution  can 
be  obtained  from  an  elementary  magneto-hydrodynamic  analysis  as 
well  as  from  an  asymptotic  Boltzmann  equation  analysis,  but  this 
question  will  be  left  to  a  later  publication. 
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5.   THE  CONTINUUM  SOLUTION 

We  consider  the  equations  of  motion  of  a  compressible  gas 
linearized  about  a  constant  density  p^,  pressure  p^,  and  velocity 
u.  : 


r  , 


i^n'    V)p  +  Pn  div  U  =  ap  U 


o  o 


(3.1)    <  Po(Uq-  V)u  +  Vp  =  -Uo(^Po%) 


(u^.  V)(P-ao  P)  =  'YPo^r 
V. 

For  generality,  arbitrary  sources  of  mass,  momentum,  and  entropy 
have  been  included;  they  will  be  specified  later.   The  momentxim 
source  is  taken  in  the  direction  of  the  unperturbed  flow.   The 
quantity  a^  is  the  sound  speed  of  the  unperturbed  flow.   The 
normalization  factors,  p^  and  u^,  have  been  inserted  for  simplicity 
of  the  ensuing  formulas;  a,  3,  and  7  are  dimenslonally  reciprocal 
lengths. 

An  equation  for  p  alone  is  obtained  by  subtracting  the  mass 
and  momentum  equations  after  operating  with  (u^-  V)  and  div 


respectively;  p  is  eliminated  by  use  of 


(3.2) 
where 
(3.5) 


P  =  ^o  P  +  Po^o'Y* 


■/ 


7dx 


-00 


(we  take  the  x-coordinate  in  the  direction  u^).   This  equation  is 

'^O 

where 
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A  = 


+ 


+ 


Sx' 


Sy^ 


5z' 


(3.5){n  =  (1-M^)  "^ 


$x 


2 


+ 


Sy^ 


V 


o   o 


S' 


Sz' 


It  Is  convenient  to  write 

(5.6)  Ay*  =  [^7*  +  M^  ^y/?^x 
In  i^A)   to  obtain 

(3.7)  D(p/Po  +  ^^'>'*)  "  -l^^  ^l^^  +  P  +  '^^'V) 
We  Introduce  a  solution  J2f  of  the  auxiliary  equation 

(5.8)  1^0   =  a  +  P  +  M^Y 
In  terms  of  which  p  is  given  by 

(3.9)  p/Po  =  -m2(7*  +  ||) 
The  velocity  is  then  found  from 

Po^o  |^=  -VP  -  ^o(^Po^o) 

2       2      -^ 
=  -  V(a  p  +  P^u^Y*)  -  u  (6p  u  ) 
v^   o^   '^o  o       o^  ^o  o' 


or 
(5.10) 

where 
(5.11) 


=  p  u^  ^  -  u  (pp  u  ) 
^o  o  dx    o^^^o  o' 


u  =  u^VJ^  -  UqP* 


&' 


X 

=  I  3dx, 
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The  relevant  singularities  are  6-functlons  for  the  momentvim  source 
and  the  entropy  source 


(3.12) 


p  =  I3^6(x,y,z) 
y   =   73_5(x,y,z) 


3-,  and  y-,  are  dimenslonally  areas. 

For  the  mass  source,  the  singularity  might  be  a  dlpole 


(3.15) 


a  =  a 


1  "^ 


5(x,y,z)  , 


but  we  shall  omit  this  for  the  present. 

We  must  now  distinguish  subsonic  flow,  M  <  1,  from  supersonic 
flow,  M  >  1.   In  the  former  case,  M  <  1,  we  introduce  the  fundamental 
solution  of  the  distorted  Laplace  equation 

(3.1M  i    i  =   -  T^ 

r^  =  x^  +  {l-M~){y^+z^). 

With  the  combined  momentum  and  entropy  source,  we  take 
0  =   (3  +  M  7,  )-i^  and  the  solution  is 

(5.15)  ^  P  -  -Po^o^^l  ■^  M^'Yi)c^l^^x 

u  =  -PiU^6*  +  u^(3^  +  M^'Yi)V^ 
The  expression 

X 

(5.16)  6*  =  r  6(x,y,z)dx  =  H(x)6(y,z) 
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FIG.  la 


FIG. lb 


FIG. 2a 


FIG.  2b 
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FIG.  30 


FIG.  3b 


FIG.  40 


FIG.  4b 
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represents  a  singularity  on  the  positive  half  of  the  x-axls, 
X  >  0.   It  Is  a  5-functlon  in  y  and  z  for  x  >  0  and  It  Is  zero 
for  X  <  0;  (H(x)  Is  the  Heavlslde  function).   This  singularity 
represents  a  wake  behind  the  object.   The  entropy  source  produces 
a  wake  In  the  density  (also  the  temperature)  and  the  momentum 
source  produces  a  wake  In  the  velocity.   The  pressure  Is  contin- 
uous (except,  of  course,  at  the  origin).   The  regular  part  of  the 
solution,  expressed  In  terms  of  i/ ,    Is  analytically  the  same  for 
the  momentiim  source  and  for  the  entropy  source.   In  Fig.  la  we 
show  the  contours  r^.  =  constant.  In  Fig.  2a  the  contours  p  =  constant, 
and  In  Fig.  5a  the  stream  lines  of  u.   It  Is  Interesting  to  note 
that  the  mass  flow  vector  field,  p  u  +  pu  ,  Is  the  superposition 
of  a  radial  (Isotropic)  source  and  a  flow  on  the  positive  x-axls, 
Fig.  4a.   This  subsonic  flow  is  not  symmetric  as  are  many  solutions 
of  invlscid  problems;  the  reason  is,  of  course,  the  momentum  source 
which  distinguishes  upstream  from  downstream. 

For  the  supersonic  case,  M  >  1,  the  analysis  Is  quite 
similar.   We  use  the  fundamental  solution 

(3.17)     I        i'  =   -      ^ 


2^M 


=  x'^  -  (M^-i)(y^+z^) 


2 
To  be  precise,  ilf   has  the  given  value  in  the  cone  r^  >  0,  x  >  0, 

and  it  is  zero  elsewhere.   No  confusion  should  arise  from  the  use 

of  the  same  symbol  V  as  before.   Taking  an  entropy  and  momentiom 

source  as  before,  we  have  0  =   {^-,    +  M  y-,)if ,   and  we  see  that  the 

solution  is  formally  identical  to  the  previous  one,  equation  (3.15)j 
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noting  only  the  change  In  the  meaning  of  if.      There  is  one  point 
of  caution  however.   This  solution  must  be  interpreted  as  a 
distribution  rather  than  as  a  conventional  function.   It  is  not 
only  the  appearance  of  the  term  6*,  but  also  the  non-integrable 
expressions  J^t/'/^x  and  yV^  which  must  be  interpreted  as  distribu- 
tions.  Moreover,  the  singularities  of  these  functions  (which 
occur  on  the  characteristic  cone,  rj^  =  0)  are  not  of  the  simple 
type  that  can  be  expressed  as  a  6-function  or  its  derivatives. 
These  complications  disappear  if  the  source  terms  are  smeared 
out  to  a  finite  size.   One  must  remember  to  interpret  the  solution 
in  terms  of  distributions  when  performing  the  smearing. 

The  contours  rj^,  =  constant  and  p  =  constant  are  shown  in 
Figs,  lb  and  2b  respectively,  and  the  velocity  field,  u,  and  mass 
flow  field,  p  u  +  pu  ,  are  given  in  Figs.  3b  and  4b  respectively. 
The  mass  flow  field  is  radial,  but  it  is  not  isotropic  as  it  was 
In  the  subsonic  case. 

It  must  be  recalled  that  these  solutions  are  the  perturbations 
over  the  original  p  and  u  .   It  is  somewhat  disturbing  to  obtain 
a  "small"  perturbation  which  is  not  only  infinite  at  the 
characteristic  cone  but  even  non-integrably  infinite,  but  this 
difficulty  is  resolved  when  we  consider  the  finite  size  of  the 
source. 

4.   SOLUTION  USING  THE  BOLTZMANN  EQUATION 

The  physical  problem  is  the  flow  around  an  object  which  is 
small  compared  to  the  mean-free-path.   Since  the  latter  is  the 
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only  length  parameter  which  appears  in  the  Boltzmann  equation, 
we  can  replace  this  object  by  a  point  source.   We  take  a  Boltzmann 
equation  with  modified  collision  term  (Ref.  (11))  and  a  source 
term 


(4.1)      e-  H  =  v(f^  -  f)  +  a-(06(x) 


The  collision  frequency,  v,  is  taken  to  be  a  constant.   The 
distribution  function  f(|,x)  is  normalized  as  a  mass  density. 
We  have 

p  =  J  fd^ 

(4.2)  ^  p  u  -  J  ^fd^ 
p  =  pRT  =  J  |(e   -  u)^  fde. 

We   Introduce   the   notation 

(4.3)  f^°)(0    =  ^°     V?  exp[-(^-uj2/2RT^] 

(27rRT    )^'^  o  o 

^  o 

for  the  unperturbed  distribution  at  infinity  and 

(4.4)  f_(e,x)  =  P   v'2  exp[-(^-  u)V2RT  ] 

°        (27rRT)^/^  ° 

for  the  locally  Maxwellian  distribution  with  parameters  p",  u,  f 
obtained  from  f(^,x)  at  the  given  point  x.   The  source  cr(^)  is 
assumed  to  be  given  and  is  easily  related  to  the  given  object. 
Specifically,  the  radial  flow  which  arises  from  the  source  cr(4) 
is  identified  with  the  free  flow  perturbation  solution;  the  cone 
on  which  this  perturbation  differs  from  zero  has  shrunk  to  a  ray. 
In  the  case  of  a  body  which  is  symmetric  about  the  axis  u  ,  we 
have  no  lift  force  and  we  can  take  the  same  symmetry  for  the 
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source  o-  (^ ) 


.(o) 


We  linearize  about  the  absolute  Maxwellian  f^  '   as  follows 

r 


.(o)  ^ 


(4.5)   1 


P  =  Po  +  P 


u  =  u^  +  u 
P  =  Po  +  P 

T  =  T   +  T 
o 


Introducing  the  relative  velocity 


(4.6) 

c   =   ^    -  u^ 

we  have 

P   =  J    gd| 

r 

PqU  +  pu^   =  J    ^gd^ 

(4.7)       \ 

p^u  =  1  cgd^ 

f  1     2 
P  =  J   5  c  gd« 

p^RT  =    r  (i  c^   -   RT^)gde 

Consistent  with  the  linearization,  we  expand  f  , 


(4.8) 


=  f<°'    1  +  ^.  (.S^-I)S  -^ 


?RT~   2'  T^  "  RT^ 
o       o     o 


and  insert  this  in  (4.1) 


(4.9)   e-||.vg  =  vf(°)|£-.(5^-i)|-.|T;j 


+  cr6. 


The  significance  of  this  linearization  requires  some  elabora- 
tion.  It  was  stated  previously  that  the  free  flow  perturbation 
solution  decays  as  l/r^  from  the  object.   This  is  true  of  any  of 
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the  macroscopic  quantities,  p,  u,  T,  but  not  of  f  itself.   For 
f  we  have  a  perturbation  of  fixed  magnitude,  but  one  which  is 
confined  to  a  smaller  cone  as  r  Increases.   However,  we  point 
out  that  the  particular  form  (4.1)  of  the  Boltzmann  equation  is 
linear  in  f;  the  nonlinearity  occurs  only  in  the  appearance  of 
the  macroscopic  quantities  p,  u,  T  in  the  local  Maxwellian  f  . 
Hence  the  linearization  applies  beyond  several  object  diameters. 
Next  we  Introduce  the  (largely  dimenslonless) quantities 


r  ■: 


p  =  p/Pq,   T  =  T/T^,   p  =  p/p^ 

u  =  u//w- ,   u^  =  u^  //rt: 


e  =  ^//^  ,   c  =  c//RT^  =  ^  -  ^o 


(4.10)  \ 


(RTJ^^^  g/p 


o 


OD  = 


^72 


exp 


(-  \c^)    =  (RT  )5/2  ^(o)/p 


o 


(27r)- 

^  =  RTq  ct/p^ 

1/L  =  v//RT^  . 

The  only  dimensional  quantities  remaining  are  o-  (a  reciprocal 
length)  and  the  mean-free-path,  L.   Equation  (4.9)  takes  the  form 


(4.11) 


4-   §e  +  T-  g  =  r  IP  +   (2  ^     -  f)T  +  c-u  I 


^x    ■    L  ^       L 
Instead  of   (4.7)   we  have 

r 


+    0-5 


(4.12)        ) 


P   =  I  gdl 
u  =    I    cgdC 

=  J  i  c^gde 

=  j   (i  c2   -   l)gd4 
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We  solve  the  integro-dlfferential  equation  (4.11)  by  Fourier 
transform;  without  confusion  we  may  use  the  same  symbol  for  the 
function  of  x  and  its  transform, 

(4.13)     0{x)    ->    (   e^^'^  0(k)dk. 

For  the  Fourier  transformed  equation,   we   obtain 


(4.14)        g(|,k)    =  -^i^  |p   +   (|  S^   -  |)T  .   c.u]  +  - 

1+iU-k    I  J       1 


+iU-k 


The  Fourier-transformed  distribution,  g,  is  explicit  as  a  function 
of  velocity,  and  as  a  function  of  k  it  is  a  linear  combination  of 
its  moments  p,  T,  and  u  plus  a  source  term.   If  we  multiply  this 
equation  by  1,  c,  and  (^  c   -  1)  successively  and  then  integrate 
over  velocity,  we  obtain  a  set  of  five  linear  inhomogeneous  equa- 
tions to  solve  for  p,  T,  and  u.   The  ultimate  solution  will  then 
be  obtained  by  inversion  of  the  Fourier  transforms. 
The  five  linear  equations  are 

--/  fidi        -  I  |(c2-3)fidf   -|  CgfidI       \ 

(4.15)  [  -|  5(c^-5)M?  l-j"  ^(c2-3)2fid|'  -||c^(2r2_3)fid| 

-J  c^fid?       -   j  ^^{c^-3)mi      1-j  c^^gfidf   /  \u 


=   L 


where 
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r 


Q.  = 


CD 


■^       J    l+iL?- 


(4.16)     \ 


'^2   = 


Ll-k 
l+lU-k 


(T-z 


_  r  c  §^dT 


V 


^        J    l+iL^-k 

The  matrix  in  (4.15)  has  five  rows  and  columns  but  is  written 
in  an  abbreviated  notation  in  which  the  third  "vector"  row  and 
column  represent  three  rows  and  columns  respectively  with  a  3x3 
block  as  the  lower  right  entry.   To  perform  the  Fourier  inversion 
explicitly,  we  consider  the  limit  of  small  L  or  large  x;  these 
are  equivalent  since  all  flow  quantities  are  functions  of  x/l. 
To  lowest  order,  approximating  n  by  co,  the  determinant  of  the 
system  (4.15)  is  found  to  be  zero.   This  is  consistent  with  the 
fact  that  the  inhomogeneous  terms  ct; ,  a^,   S%  are  multiplied  by 
the  small  quantity  L.   To  the  next  order  we  write 

(4.17)  n  =  aj(l-iLl'k)  =  co(l-iLc-k  -  iLu  -k). 

It  is  a  simple  matter  to  verify  that  this  approximation  to  fi 

yields  the  correct  asymptotic  approximation  to  the  various  integrals 

that  must  be  taken  over  n.   Inserting  this  approximation,  the 

,^      ^      A/ 

linear  system  for  p,  T,  u  becomes 
/iu^-k     0        ik 
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This  matrix  is  easily  inverted,  and  we  find 


3   ^  o   ^ 


-k' 


(4.19)   D   T 


rJ 


2,2     1  2  /-.   ,  ^2 
-  ^  k     k  -(u^-k) 


(u^-k)k^    (u^-K)l<r 


(u  •k)k^ 
^  o    s 


(a„-k)k3 


[|k2-(u„-k)2l5 


o 

-s-  k 
3  r  s 


rs 


-  5k  k 


Here  D  is  the  determinant  of  the  matrix  in  (4.l8), 

(4.20)    D  =  i(uQ-k)[|  k^  -  (a^-k)2]. 

v;e  notice  that  D  represents  the  Fourier  transform  of  the 
product  of  the  operator  Q  with  (u  •  V)-   This  operator  is  easily 
inverted  without  recourse  to  Fourier  transformation,  so  it  remains 
only  to  invert  the  right  side  of  equation  (4.19).   To  do  this,  we 
first  study  the  Integral 
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(^.21) 
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dk 
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This  is  the  Fourier  inversion  of  a  typical  expression  which  arises 
in  the  source  terms  a;,  CT"p,  and  ct^,  but  we  have  chosen  to  invert 
the  transform  before  integrating  out  the  velocity.   To  evaluate 
this  expression,  we  choose  the  fixed  vector  ^  as  a  coordinate 
axis,  i   =  ii,    0,  0)  and  write  k  =  (k-j^,  k^,  k^),  x  =  (x^,  x^,  x^). 
Integration  with  respect  to  k,  yields,  for  x,  >  0 


(4.22) 


S  =  -^  e 
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dk„dk, 
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and  zero  for  x-,  <  0.   We  recognize  the  remaining  integral  as  a 
6-function  with  respect  to  Xp  and  x^  and  conclude  that 
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1   -x./L^ 
(4.23)    S  =  -i  e  ^  H(x,)5(x2,x-,)c^(e) 

Furthermore,  for  small  L,  the  first  factor  approaches  a  6-function 
in  X,.   This  result  has  a  simple  intuitive  interpretation  of  S(|,x) 
as  a  radial  flow  arising  from  the  source  function  cr(0.   This 
flow  is  exactly  the  free-flow  solution  for  the  body  which  corresponds 
to  the  source  <t  modified  by  an  exponential  decay  representing 
collisions  with  a  fixed  background  gas  of  mean-free-patn  L. 

The  various  integrals  ai,   a%,   and  at  are  obtained  from  S  by 
integration  over  velocity  after  multiplication  by  a  suitable 
polynomial  function  of  i.      The  result  is  an  approximate  5-function 
in  X  (width  on  the  order  of  L)  with  a  numerical  factor  depending 
on  the  polynomial.   First  we  remark  that,  in  order  to  conserve 
mass,  we  must  have 


(4.24)     j  a{X)dX   =  0. 


From  this  we  conclude  that  (3",  vanishes  in  the  limit  of  small  L; 
i.e.  there  is  no  B-function  left  in  this  limit.   The  Fourier 
inversions  of  ctL  and  cr.,  are  given  by 


^2 


=  I  (^  c2-i)  S(i:,x)de   -   T^  u^bix) 

(^.25) 

6^3  =  [c  S(1,x)d?   ^  T^u^6(x) 

The  vector  o^  has  the  direction  u  in  virtue  of  the  assumed  symmetry 
of  ?(0-   If  we  assume  that  the  total  energy  of  the  gas  is  conserved 
(i.e.  the  object  has  reached  an  equilibrium  temperature  and  there 
is  no  radiation)  then 
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(4.26)    (l-crii)      =     0. 
In  this  case,  since  c  =  ^ 


2£-u  +  u   ,  we  can  write 
o    o 


"^2  = 


=  -||(^"o) 
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but  also 


from  which  we  conclude  that 


(4.27) 


Tg  -  -  3  -^3. 


Returning  to  (4.19)  and  adopting  the  6-function  approximations 
to  &"-.,   (Tp,  and  at,  we  obtain,  after  inverting  the  Fourier  transforms, 

(4.28)  ^(u^-V)nT   =  [M^T^  □  -  f(M^2  -  ^3)(^o'  "^^^^^ 


(u^-V)Qu   =  [t^u^  D  +  (M  ^2  -  T^)(v^^.  V)V]5 
We  have  used  the  relation 
(4.29)     ^0=%/^^  =  /573  M 

Solving  these  differential  equations  Just  as  in  the  previous 
section,  we  obtain  the  identical  results  that  were  obtained  there 
(viz.  equation  (3.15))  provided  that  we  make  the  identifications 


(4.30) 


'^1  =  -  -^3/^0 


7i  =  Tg/u;. 
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The  restriction  (4.27)  which  arises  from  energy  conservation  relates 
the  momentum  and  entropy  sources. 
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This  relation  is  equivalent  to  the  thermodynamic  statement  that 
the  entropy  source  in  the  wake  is  given  by  the  irreversible  work 
corresponding  to  the  product  of  the  drag  by  the  free  stream 
velocity,  u^. 

It  is  also  easy  to  demonstrate  another  universal  relation, 
viz.  that  the  dipole  component  of  the  mass  source  is  simply 
related  to  the  drag,  whatever  form  aii)   may  take. 

5.   SUMMARY  AND  CONCLUSIONS 

In  the  asymptotic  approximation  for  small  mean-free-path 
which  we  found  to  yield  the  continuum  result,  the  finite  size  of 
the  mean-free-path  enters  in  two  distinct  places.   First  there 
is  the  finite  extent  of  the  source  terms  cr,  ,  a^,  a^.   Then  there 
is  the  modification  of  the  equations  of  motion  (i.e.  of  the  matrix 
of  the  linear  system  for  p,  T,  u)  to  include  dissipation. 

The  significant  fact  about  the  mass,  momentum,  and  entropy 
sources  is  that  their  extent  is  measured  by  the  mean-free-path 
rather  than  by  the  object  size  (which  is  even  zero  in  this 
formulation).   The  sources  c,,  a^,   a--,   are  finite  in  extent  even 
though  the  original  source,  a  6,  in  the  Boltzmann  equation  is  not. 
Both  the  shock  and  the  wake  are  opened  up  to  at  least  a  mean-free- 
path  in  thickness.   Also,  the  structure  of  these  sources  yields 
a  definite  internal  structure  for  the  shock  and  for  the  wake.   This 
refinement  can  be  easily  computed  in  terms  of  the  fundamental  solu- 
tion 1^  (equations  (3.14)  and  (3.1?))  of  the  non-dissipative  theory. 
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This  elaboration  of  the  elementary  results  will  yield  a  mass 
source  which  is  different  from  zero;  of  course  the  integrated 
source  remains  zero,  but  it  is  only  in  the  limit  L  =  0  that 
there  is  no  mass  source  at  all. 

Incidentally,  if  the  magnitude  of  the  source,  o~,  is  related 
to  an  assumed  object  of  size  a  which  is  small  compared  to  L, 
then  the  actual  disturbance  in  the  shock  wave  and  in  the  wake 
will  be  small  (at  most  of  order  a  /L  )  once  these  boundary  layers 
are  smeared  out  to  the  thickness  L. 

The  finite  extent  of  the  sources  yields  a  minimum  thickness 
for  the  shock  and  for  the  wake.   In  addition,  one  naturally 
expects  that  keeping  a  finite  mean-free-path  in  the  rest  of  the 
theory  will  produce  a  progressive  widening  of  these  boundary 
layers.   Once  the  boundary  layers  thicken  to  more  than  a  few 
mean-free-paths,  it  is  natural  to  expect  the  usual  dissipative 
continu\am  theories  to  apply.   As  a  matter  of  fact,  it  is  easy 
to  verify  that  the  stresses  and  heat  flow  computed  for  the  zero 
order  distribution  function  found  in  the  last  section  are  exactly 
the  Navier-Stokes  values.   Of  course,  the  exact  nature  of  the 
"next"  approximation  depends  on  how  one  decides  to  improve  the 
Fourier  inversion.   Nevertheless,  one  can  confidently  expect  these 
boundary  layers  to  widen  parabolically  as  would  be  predicted  by 
the  usual  viscous  and  heat  conduction  laws.   This  widening  is  the 
ultimate  Justification  for  the  non-dissipative  solution  as  a 
5-function  in  the  wake  and  a  sharp  singularity  on  the  Mach  cone. 
On  a  large  scale,  a  parabolically  widening  boundary  layer  will 
appear  as  a  singular  line  or  surface.   To  be  precise,  at  a  fixed 
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position  in  space,  the  dissipative  solution  will  approach  the 
non-disslpatlve  one  as  L  ->  0  (the  exact  solution  is  a  function 
of  x/l).   The  strength  of  the  singularity  6*  in  the  wake  is  a 
measure  of  the  integrated  disturbance  over  the  entire  wake.   The 
flow  near  the  regions  of  nonuniform  convergence  can  be  found  by 
allowing  (x,y,z)  to  approach  infinity  on  an  appropriate  parabolic 
path  instead  of  radially;  this  has  not  yet  been  done. 

It  should  be  remarked  that  the  singularity  at  the  Mach  cone 
should  not  really  be  called  a  shock.   By  the  term  "shock"  one 
usually  refers  to  a  wave  which  has  reached  an  equilibrium  against 
the  competing  effects  of  dissipation  tending  to  widen  it  and 
nonlinearities  tending  to  steepen  it.   It  is  easy  to  verify  that 
the  strength  of  a  wave  lying  on  the  Mach  cone  decays  at  least  as 
fast  as  l/r.   This  would  imply  a  shock  thickness  increasing  at 
least  as  r,  i.e.  conically.   But  the  maximum  thickness  caused  by 
dissipative  broadening  is  parabolic  and  has  the  order  /r  .   Thus 
the  broadening  process  is  too  slow  to  ever  develop  into  a  true 
shock.   Alternatively,  one  can  say  that  nonlinearity  becomes 
increasingly  negligible  with  distance  from  the  origin.   Also, 
the  terminology  "detached"  shock  is  used  only  because  of  a 
general  qualitative  similarity  between  the  flows.   For  a  true 
shock,  the  strength  directly  ahead  of  the  object  is  determined 
by  the  Mach  number  alone.   In  the  flows  we  describe,  with  the 
object  size,  a,  small  compared  to  L,  the  strength  of  the  distur- 
bance ahead  of  the  object  cannot  greatly  exceed  the  small  quantity 
a  /L  .   Nevertheless,  one  can  expect  a  continuous  transition  of  the 
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flow  described  here  into  a  true  continuum  flow  with  a  detached 
shock  as  a/L  increases. 

We  return  to  *"-he  question,  how  to  improve  the  elementary 
results  given  above.   In  the  neighborhood  of  the  object,  the 
conventional  free  flow  solution  is  accurate.   Far  away  and  outside 
the  boundary  layers  the  elementary  continuum  flow  is  valid.   One 
point  is  abundantly  clear.   The  direct  continu\im  analysis  is 
much  simpler  than  the  equivalent  asymptotic  limit,  L  ->  0  even 
with  the  simplified  Boltzmann  collision  term  which  was  adopted. 
One  might  wish  to  improve  the  continuum  model  by  taking  the 
Navier-Stokes  or  thirteen  moment  equations,  judiciously  modify- 
ing the  source  terms  as  suggested  by  the  Boltzmann  equation  analysis. 
Several  steps  seem  natural.   The  first  is,  as  mentioned  above,  to 
use  the  same  non-dissipative  equations  but  with  the  smeared  out 
source  terms;  this  answer  is  explicit  as  an  integral  over  the 
source  term  with  the  fundamental  solution,  ilf ,   as  a  kernel.   Next 
one  might  try  the  Navier-Stokes  equations  with  the  same  source  terms; 
or  the  thirteen  moment  equations,  incorporating  additional  sources 
in  the  stress  and  heat  equations  by  taking  second  and  third 
moments  of  the  function  S  (equation  (4.25))-   Either  of  these 
can  be  expected  to  give  results  which  are  reasonably  accurate 
over  the  entire  range  of  distances. 

There  is  one  other  type  of  improvement  which  we  suggest   for 
its  inherent  illumination  of  the  problem  rather  than  for  its 
practicability.   The  same  function  S  which  arose  in  the  analysis 
of  the  asymptotic  behavior  as  L  ->  0  arises  in  another  context, 
viz.  in  the  treatment  of  the  free  flow  problem  near  the  body,  which 
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can  be  considered  as  the  asymptotic  behavior  as  L  ->  oo  I   The 
function  S  Is  exactly  the  free  flow  perturbation  (i.e.  the 
difference  between  the  unperturbed  flow  and  the  free  flow  solutlon- 
thls  Includes  the  removed  particles  in  the  shadow  as  well  as  those 
reflected  by  the  body)  modified  by  colllslonal  losses  with  the 
oncoming  stream.   This  function  can  be  interpreted  as  the  first 
step  in  an  iteration  process;  the  next  step  would  include  as 
sources  the  particles  whose  loss  by  collision  is  Inherent  in  the 
exponential  decay  of  the  function  S.   One  could  then  choose  a 
higher  order  step  in  this  Iteration  about  the  free  flow  solution 
as  an  appropriate  source  term  instead  of  S  itself.   In  this  way, 
the  solution  would  be  simultaneously  improved  from  both  ends,  near 
the  body  and  far  from  it. 

This  procedure  can  be  generalized.   Let  us  suppose  that  the 
distribution  function  is  known  exactly  (or  to  some  approximation) 
in  any  region,  small  or  large,  which  surrounds  a  source  or  a 
finite  body.   We  replace  this  region  by  a  fictitious  distributed 
source  located  on  the  boundary  of  the  region  such  that  the  source 
exactly  duplicates  the  molecular  inflow  and  outflow.   At  a  large 
distance  from  this  region,  the  source  can  be  considered  to  be  a 
point  source  in  the  Boltzmann  equation,  i.e.  it  is  exactly  the 
function  cr  .   It  is  clear  that  the  same  argument  applies  to  the 
replacement  of  a  rigid  object  of  any  size  by  an  appropriate  source. 
The  difficulty  is  that,  if  the  object  is  not  in  free  flow,  the 
equivalent  source  term,  cr,  is  not  known  but  depends  on  a  solution 
for  the  distribution  function  near  the  object.   To  the  lowest 
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order,  however,  far  from  the  body,  we  have  seen  that  a  single 
number,  the  value  of  the  momentum  source  or  drag,  suffices  to 
specify  the  asymptotic  flow.   V/e  conclude,  therefore,  that  the 
elementary  solution  obtained  in  Section  3  represents  the  correct 
asymptotic  flow  far  from  an  arbitrary  body  of  any  size  and  speed. 

Finally,  we  conclude  with  the  claim  that  an  entirely 
equivalent  analysis  including  all  of  the  complexities  associated 
with  an  ionized  conducting  gas  seems  to  be  quite  feasible.   In 
particular,  one  can  expect  that  the  simplest  non-dlssipative 
continuum  equations  will  give  the  correct  asymptotic  behavior. 
Unfortunately,  this  continuum  theory  can  only  be  solved  explicitly 
in  certain  special  cases,  but  approximation  techniques  to  cover 
the  entire  problem  are  being  developed  and  should  be  at  hand 
shortly. 
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